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Abstract 

A combination of algebraical and numerical techniques for calculating two-loop 
top-quark and Higgs-boson corrections to electroweak precision observables 
like Ar or the p-parameter is presented. The renormalization is performed 
within the on-shell scheme. The results of the calculations are valid for arbi- 
trary values of m t , M H and of the gauge-boson masses. An example is treated 
where the full result is compared to the result obtained via an expansion up 
to next-to-leading order in m%. As an application, results for the Higgs-mass 
dependent top-contributions to Ar are given. 
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A combination of algebraical and numerical techniques for calculat- 
ing two-loop top-quark and Higgs-boson corrections to electroweak 
precision observables like Ar or the p-parameter is presented. The 
renormalization is performed within the on-shell scheme. The re- 
sults of the calculations are valid for arbitrary values of nit, -Wh 
and of the gauge-boson masses. An example is treated where the 
full result is compared to the result obtained via an expansion up 
to next-to- leading order in m t . As an application, results for the 
Higgs-mass dependent top-contributions to Ar are given. 

1 Introduction 

After the discovery of the top-quark the Higgs-boson remains the only miss- 
ing ingredient of the (minimal) electroweak Standard Model (SM). At the 
moment, the mass of the Higgs-boson, Mh, can only very mildly be con- 
strained by confronting the SM with precision data. In order to improve on 
this situation, a further reduction of the experimental and theoretical errors 
is necessary. Concerning the reduction of the theoretical error due to miss- 
ing higher order corrections, in particular two-loop top-quark corrections and 
higher order QCD corrections to the quantity Ar derived from muon decay and 
to the p-parameter are of interest. Recently the leading two- loop top-quark 
and Higgs-boson corrections to these quantities [1] have been supplemented 
by the calculation of the full Higgs-boson dependence of the leading m\ con- 
tribution [2] and furthermore by inclusion of the next-to-leading top-quark 
contributions [3]. Since both the Higgs-mass dependence of the leading m\ 
contribution and the inclusion of the next-to-leading term in the m t expan- 
sion turned out to yield sizable contributions, a more complete calculation of 
the top-quark corrections would be desirable, where no expansion in m t or 
M H is made. 
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Fig. 1. One-particle irreducible topologies of two-loop self-energies which do not 
factorize into one- loop contributions. 

In this article we describe techniques with which such a calculation can be 
carried out. The calculation of top-quark contributions to Ar and other pro- 
cesses with light external fermions at low energies requires in particular the 
evaluation of two-loop self-energies at non-zero external momentum, while 
vertex and box contributions can mostly be reduced to vacuum integrals. We 
perform the renormalization within the complete on-shell scheme, i.e. we use 
physical parameters throughout. All calculations are performed in an arbitrary 
i?£ gauge with free gauge parameters i — A,Z,W '. 

2 Algebraic Evaluation 

The problem of evaluating top-quark or Higgs-boson corrections to four-fermion 
processes with light external particles at low energies in the on-shell renormal- 
ization scheme requires the calculation of two-loop self-energies with massive 
particles at non-zero external momentum. The algebraic evaluation of these 
contributions can be carried out in a highly automatized way by means of the 
Mathematica packages FeynArts 2.0 [4] and TwoCalc [5]. 

The package FeynArts 2.0 is used for generating the relevant one-loop and 
two-loop diagrams and counterterm graphs. For this purpose, we have imple- 
mented the complete one-loop counterterm contributions of the SM (using the 
complete on-shell renormalization scheme) as well as the appropriate two-loop 
terms into the model file of the package. 

The package TwoCalc reduces general two-loop self-energies to a minimal set 
of standard scalar two-loop integrals of the form 

T il i 2 ...i l (p 2 ;rnl,ml, . . . ,mf) = 

d D q 1 f d D q 2 1 



f d gi f 
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« 2 (2^)°- 4 J in 2 (2TTLi) D - 4 [kl - m\] [kf 2 - m|] • ■ • [fc? - mf 

which correspond to the basic two- loop topologies depicted in Fig. 1. In (1) 
ki t denotes the momentum of the i-th propagator and mi its mass. The k^ are 
related to the integration momenta q\ and q 2 and the external momentum p 
via ki — qi, k 2 = qi + p, fc 3 = q 2 — qi, k 4 = q 2 , k 5 = q 2 + p. The topologies 
shown in Fig. 1 correspond to the scalar integrals T12345, T11234, ^1234, ^234, and 
Tiiu, respectively. The analytical expression for T11234 can be obtained from 
Tus4 by partial fractioning or taking the derivative with respect to m\. Other 
integrals with higher powers of propagators are treated in the same way. For 
the general case one therefore needs to evaluate only four different types of 
two-loop scalar integrals. 

The algebraic evaluation with TwoCalc is performed for arbitrary values of all 
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particle masses, the invariant momentum p 2 and the space-time dimension D. 
A general gauge is used which is specified by one gauge parameter £j, % — 
A, Z, W, for each vector boson. The algorithm of TwoCalc has been described 
in detail in Refs. [5]. For the contraction of Lorentz indices, reduction of the 
Dirac algebra and evaluation of Dirac traces, which can be worked out like 
in the one-loop case, the routines of the program FeynCalc [6] are used. The 
further evaluation is based on a method for the tensor integral decomposition 
of two-loop self-energies [5] and furthermore makes use of certain symmetry 
properties of the two-loop integrals. 

The counterterm contributions are also calculated with TwoCalc. The sum of 
unrenormalized two-loop diagrams and diagrams with counterterm insertions 
is expressed in terms of a minimal basis of standard scalar integrals. This is 
a very convenient feature, since it allows to directly check at the algebraic 
level whether the gauge-parameter dependence of the result drops out. As 
a further nontrivial check which can directly be read off from the algebraic 
result, all field renormalization constants of internal particles must drop out in 
the sum of all contributing diagrams. After inserting the divergent part of the 
two-loop integrals and decomposing the one-loop integrals into divergent and 
finite parts, also the UV-finiteness of the result can be checked algebraically, 
i.e. all terms proportional to 1/5 2 and 1/5, where 5 = (4 — D)/2, must cancel. 

3 Two-loop On-shell Renormalization 

For our calculations we use the complete on-shell renormalization scheme at 
the two-loop level, i.e. we use as parameters the masses of the physical particles 
and the electric charge e, and we introduce field renormalizations such that the 
residues of all propagators are equal to one (see e.g. Ref. [7]). The formulation 
in terms of physical parameters is particularly important in view of calculating 
corrections to the quantity Ar derived from muon decay, which relates the 
Fermi constant to the W-boson and Z-boson masses, M\y and Mz, and to 
the electromagnetic fine structure constant a = e 2 /(47r). 

We perform the renormalization in such a way that the gauge-fixing term in 
the Lagrangian does not give rise to counterterm contributions, i.e. the renor- 
malization of the parameters and fields in the gauge-fixing term is canceled 
by appropriate renormalizations of the gauge parameters. It should be noted 
that the renormalization prescription chosen for the gauge-fixing term deter- 
mines the renormalization of the ghost sector, which at the two-loop level also 
enters the quantities in the physical sector of the theory. According to the 
described renormalization procedure we have derived the counterterms for the 
two-point and three-point vertices containing ghost fields. For the example of 
the complete two-loop top-quark contribution to the W-boson self-energy we 
have checked that insertion of these renormalization constants does in fact 
yield a finite result for the two-loop self-energy. 

As a specific example we briefly consider the Higgs-dependent top-quark cor- 
rections to Ar. It is easy to see that Higgs-dependent top-quark corrections 
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can, except for the renormalization, enter Ar only via two-point functions. Be- 
sides these self-energy corrections one also has to consider contributions aris- 
ing from the renormalization of the relevant one-loop and two-loop three-point 
functions. The two-loop renormalization constants receiving Higgs-dependent 
top-quark corrections are the charge renormalization constant, 5Z e ^ 2 ), and the 
counterterm of the electroweak mixing angle, Ss 2 ^ r 2 -\, while the field renormal- 
ization constant SZ W ^ 2 ) (and also SZ W ^) cancels in the sum of the contribut- 
ing diagrams. 

4 Decomposition and numerical Evaluation of scalar Integrals 

It is known since several years that massive two-loop integrals are in general 
not expressible in terms of polylogarithmic functions [8]. Our approach for a 
numerical evaluation of scalar self-energy integrals are one-dimensional inte- 
gral representations. They allow for a very fast calculation of these functions 
with high precision. 

Such a numerical approach can of course only be applied to finite functions, 
while UV- and IR- divergencies of the two-loop integrals have to be kept under 
control. We decompose the integrals into divergent and finite parts and check 
algebraically that the divergencies cancel. For some integrals we make use of 
the fact that the two-loop vacuum integrals can be calculated analytically 
(see e.g. Ref. [9]). Then we numerically calculate suitable combinations of the 
non-vacuum and the vacuum integrals. For some diagrams it proved useful to 
separate the two-particle cut and the three-particle cut contributions [10], or 
to combine both methods. As an example the momentum derivative of T11234 
can be calculated as 

^T 1123 4(p 2 , ml, ml, m\, mj) = ^Bo{m\, m\, m 2 4 )-^B (p 2 , m\, m%)j 

+DT 112UC 3(P 2 , ml, m 2 2 , m\, m 2 4 ) , (2) 

where DT n234C 3, the three-particle cut contribution, is finite, while the di- 
vergencies are contained in one-loop I?o-integrals, which can be calculated 
analytically. 

For the numerical evaluation there remains an expression for the finite part 
which is in general very extensive. The main ingredients of the numerical 
evaluation are one-dimensional integral representations which we have imple- 
mented in C++. For the integration we apply an adaptive Gauss-Kronrod 
algorithm [11]. 

For those topologies which contain a one-loop self-energy insertion we make use 
of a dispersion representation of the subloop which has been discussed in detail 
in Refs. [10]. In these representations the two-particle cut contributions have 
been separated as products of i?o-functions. All momentum or mass derivatives 
of the functions T 234 and T 1234 can be calculated starting from these formulae. 
For the master topology, X12345, we refer to the formulae presented in Ref. [12]. 
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The basic functions in our integration kernels are .Bo-functions and the discon- 
tinuities of B$-{ unctions. In some cases it is useful for the numerical stability to 
subtract the asymptotic behavior of the 5 -function for a large mass- variable. 
It is given by 

1 S 7T1 2 2 2 

B (p 2 ; s, m 2 ) = - + 1 - log + — log — + — + B 0rest (p 2 ; s, m 2 ) . (3) 

Calculating physical processes we encountered in some parameter regions huge 
cancellations among the contributions of the diagrams or the scalar integrals. 
Therefore we perform in these cases the calculations inside the integral repre- 
sentations with quadruple precision [13]. Typical computation times are 0.03 
seconds for the evaluation of the functions T234 and T1234 and 0.3 seconds for 
the evaluation of T12345 to ten digits precision on a workstation DEC 3000 
AXP. Using quadruple precision slows down the calculations by a factor of 
about 10, which results e.g. in 40 seconds for the calculation of the Higgs- 
dependent top-quark contributions to the W-boson self-energy for one set of 
parameters. 

5 Results 

5. 1 Comparison with the expansion in m t 

In this section we compare our results for the two diagrams shown in Fig. 2 
with the results of an expansion up to next-to-leading order in (see 
Ref. [3]), which takes into account terms of order m\ and m 2 . This expansion 
is performed in two regions, namely in the light Higgs region ("light Higgs 
expansion" ) and in the heavy Higgs region ( "heavy Higgs expansion" ) [3] . 



w ; ° J W W \ / 



H 

(a) (b) 



Fig. 2. The diagrams for which we compare our results with the expansion in m t . 
We first consider asymptotically large values of m t and check whether our full 
result and the expansion in m t agree in this region. In order to compare with 
the heavy Higgs expansion, both m t and Mr are made asymptotically large, 
while for comparison with the light Higgs expansion Mh is kept fixed. Fig. 3 
shows for diagram (a) the difference between the finite parts of our result and 
the heavy Higgs expansion, divided by m t , as a function of m t . For the self- 
energies we have here and below used units of Myyfs, which amounts with 
a(M|) = to 1.52 x 10" 6 M^. The difference between the full result and 
the expansions in m t is expected to be at most of the order m t , as it is the case 
in Fig. 3. We have checked that also for the light Higgs expansion of diagram 
(a) and for both expansions of diagram (b) one finds agreement between the 
full result and the expansions in the asymptotic region. 



1 P. Gambino has kindly provided us with his results for the expansions of these 
diagrams. 
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Fig. 3. Diagram (a): Difference between full calculation and heavy Higgs expansion, 
divided by m t (p 2 = M^). 

Next we consider the physical parameter region, i.e. we put m t = 175 GeV, 
and compare the full result with the expansions for different values of M H . 
This is shown in Fig. 4 for diagram (a) and in Fig. 5 for diagram (b). We 
find relatively good agreement for diagram (b), while for diagram (a) for most 
values of Mr the relative deviations are larger. It should be noted, however, 
that the numerical contribution of diagram (b) is very large so that in the 
final result large cancellations can be expected. 
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Fig. 4. Diagram (a): Full calculation and expansions (m t = 175 GeV, p 2 = M^). 
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Fig. 5. Diagram (b): Full calculation and expansions (m t = 175 GeV, p 2 = M^). 
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5.2 Sensitivity of the two-loop top-quark corrections to Ar on the Higgs-boson 
mass 

In this section we study the sensitivity of the two-loop top-quark corrections 
to Ar to the Higgs-boson mass by considering the quantity 

A^W^h) = Arg>(M H ) - Arg } p (M H = 60 GeV), (4) 

where Ar^F(Mn) denotes the complete two- loop top-quark contribution to 
Ar. In Fig. 6 the variation of Ar^jj 5 ubtr (MH) with the Higgs-boson mass is 
shown in the interval 60 GeV < Mh < 1 TeV for various values of m t . The 
change in Ar^(M-a) induced by varying Mh in this interval is found to be 
about 0.001. It is interesting to note that the absolute value of the correction 
is maximal just in the region of m t = 175 GeV, i.e. for the physical value of 
the top-quark mass. 



-0.0002 
-0.0004 
Ar (2) -0.0006 
-0.0008 
-0.001 
-0.0012 

100 200 300 400 500 600 700 800 9001000 -f^ 
Fig. 6. Two-loop top-quark contributions to Ar subtracted at Mh = 60 GeV. 
In Fig. 7 the corresponding one-loop correction, 

^ r (l),subtr (M H ) = Ar (1) (M H ) - Ar (1) (M H = 60 GeV) (5) 

is shown together with the combined contribution Ar( 1 ) )Subtr (M H )-|-Ar|2j ) subtr (MH). 
As one can see in the plot, the one- loop and two- loop corrections enter with 
different sign, i.e. the sensitivity of Ar to the Higgs-boson mass is lowered by 
the inclusion of the two-loop top-quark corrections. The size of the two-loop 
correction Ar^j 5 ubtr (M H ) is about 10 percent of the one- loop contribution. 

6 Conclusions 

In this article we have outlined techniques for the calculation of electroweak 
two-loop corrections associated with the top quark or the Higgs boson to preci- 
sion observables like Ar. No expansion in the particle masses is performed, i.e. 
the calculations are valid for arbitrary values of m t , Mh and the gauge-boson 
masses. The calculations are performed with the help of computer-algebra 
packages, which provide a high degree of automatization, and very efficient 
C-routines for the numerics, which are based on one-dimensional integral rep- 
resentations of the two-loop scalar integrals. For an example we have compared 
our results to those obtained via an expansion in the top-quark mass up to 
next-to-leading order. We have furthermore calculated the sensitivity of the 
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Fig. 7. One- loop and two- loop top-quark contrib. to Ar subtracted at 
M H = 60 GeV. 

two-loop top-quark corrections to Ar with respect to the Higgs-boson mass 
and compared these contributions to the one- loop result. 
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